(ttolUgp  af  SIthFral  Arta 
ffithrarg 

THE  GIFT  OF   :^Vve.~l^LC^W«>^ 


BOSTON  UUIVERSITY 
GHADUATB  SCHOOL 
Thesis 

INTRODUCTION  TO  THE  SOLUTION  OF  ORDINaHI 
LINEAR  DIFFERENTIAL  EQUATIONS  WITH 
CONSTANT  COEFFICIENTS 

TTilliam  Arthur  Lowell 
(B.S.  Bates  College,  I9IS) 
suhnitted  in  partial  fulfilment  of  the 
requirements  for  the  degree  of 
Master  of  Arts 
103^ 


TABLE  OF  COITTEKTS 


Introduction 

Article  Page 

1.  Brief  History   1 

2.  Importance  .    •   1 

3.  Metliod  of  Procedure   2 

Definitions  and  Assumptions 

k.  Differential  Equation    2 

5,  Linear  Differential  Equation    2 

Sm  Homogeneous  -  Non-xiomogeneouE  Linear 

Differential  Equation    3 

7,   Solution  of  Differential  Equation  •    3 

g.  Assumptions   ••••  3 

Eoraogeneous  Linear  Differential  Equation 

9.   General  Tlieorera   3 

10.  Constant  Coefficients,     2oot3  of  Auxiliary- 
Equation  Real  and  Distinct  ,   ,  5 

11.  Hoots  of  Auxiliary  Equation  Complex 

Imaginary  and  Distinct                                             ...  7 

12.  Hoots  of  Auxiliary  Equation  Multiple    10 

Non-homogeneous  Linear  Differential  Equations 

13*   Greneral  Theorem     16 

lU,   Symbolic  Operator   I7 

15.  Operator  Method  for  Particular  Integral    I9 


lo.  Variation  of  Parameters   22 

17»  Undetermined  Coefficients    26 

Summary  and  Conclusions    .....  3^ 

Bililiograph/   36 


INTRODUCTION 

!•  Brief  History.  -    The  study  of  Differential 
Equations  began  very  soon  after  the  discovery  of  Differential 
and  Integral  Calculus,  to  which  it  forms  a  natural  sequel. 
In  1676  Newton  solved  such  an  equation  "by  means  of  infinite 
series  but  did  not  publish  his  results  until  1693*   the  same 
year  that  Liebnlz  encountered  them  in  his  work.     Prom  then  on 
progress  was  rapid.     Such  men  as  the  Bernoulli  brothers, 
Fontaine,  Buler,  Clairaut,  Taylor  and  Lagrange  became 
interested,  and  by  the  beginning  of  the  nineteenth  century 
many  important  results  had  been  established  but  the  methods 
developed  were  artificial,  broad  comprehensive  principles 
lacking.     Puch's  Memoirs  (1866-6S)  gave  the  theory  of  Linear 
Differential  Equations  its  birth,   the  Lie  Theory  of 
Continuous  Groups  (18SU)  advanced  it  and  led  Picard  and 
Vessibt  to  develop  a  theory  which  is  analogous  to  the  Galois 
theory  for  algebraic  equations, 

2.   I^mportance  of  Linear  Differential  Equations .  - 
The  problem  of  solving  equations  of  this  type  is  constantly 
met  in  scientific  investigations,  particularly  in  those 
dealing  with  the  study  of  motions  of  va.rious  kinds.     They  are 
encountered  in  Algebra,  Geometry,  Mechanics,  Physics,  Chemistry 
and  Astronomy,     Many  problems  in  tangency,   curvature,  envelopes 
oscillations  of  mechanical  systems  and  of  electric  currents, 
bending  of  beams,  conduction  of  heat,   diffusion  of  solvents. 


2  i  0  I  - 
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velocity  of  chemical  reactions,  and  celestial  mechanics 
formulate  themselves  as  linear  differential  equations  of  one 
of  two  types:   -  those  with  constant  coefficients  which  are 
the  more  important,  or  those  with  variahle  coefficients.  In 
fact,  30  great  is  the  preponderance  of  these  equations  in 
Physics  alone  that  the  study  of  them  in  this  field  might  well 
"be  called  the  study  of  Linear  Differential  Equations,  The 
frequency  of  their  occurrence  in  any  one  of  the  above  mentioned 
"branches  of  science^suf f icient  justification  for  the  special 
attention  given  them  in  this  paper, 

3,  Method  of  Procedure,  -     The  method  of  reasoning 
used  in  developing  the  subject  is  largely  the  inductive  type 
which,  I  "believe,  is  consistent  with  the  fact  that  it  is  of 
an  introductory  nature,     A"bsolute  rigor  has  "been  sacrificed 
in  some  instances  for  the  purpose  of  simplifying  the  argument, 

DEriKITIQlTS  AITD  ASSUMPTIONS 

U,  Differential  Equations,  -     An  equation  involving 
derivatives  or  differentials  is  a  differential  equation.  If 
the  equation  involves  a  single  independent  varia"ble,   it  is  an 
ordinary  differential  equation;  otherwise,  it  is  a  partial 
differential  equation, 

5,   Linear  Differential  Equation,  -     A  differential 
equation  that  is  of  the  first  degree  in  the  dependent  variable 
and  its  derivatives  is  a  linear  differential  equation.  The 
general  form  of  the  ordinary  linear  differential  equation  is 
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Where  Xq  ,  Xj_,  X^^  and  X  are  functions  of  x  or  constants, 

6.  Konogeneous  -  Non-homogeneous  Linear  Differential 
Equation.  -    A  linear  differential  equation  in  which  the 
right-hand  memljer  is  zero  is  a  homogeneous  linear  differential 
equation.     When  the  right-hand  nemher  is  a  function  of  x,  the 
equation  is  said  to  he  non-homogeneous. 

7.  Solution  of  Differential  Equation.   -  Any 
functional  relationship  (not  containing  derivatives)  between 
the  variables  which  causes  a  differential  equation  to  reduce 
to  an  identity,  is  a  solution  or  a  particular  integral  of  the 
equation.     The  meaning  of  the  term  "solution"  in  this  theory- 
is  somewhat  broader  than  that  given  it  in  Algebra,     It  is 
extended  to  include  not  only  implicit  algebraic  forms  such  as 
y  *  oC  X  yx    +  y  ,  but  also  symbolic  forms  such  as  y  =  ^f(x)dx 
which  cannot  be  expressed  without  integral  signs, 

g.  Assumptions.  -     It  is  assumed  in  this  discussion 

of  the  n*'^  order 
that  every  linear  differential  equation^has  a  general  solution, 

and  that  this  solution  involves  n  arbitrary  independent 

constants. 


HOMOOENBOirS  LIITEaR  DIFFEEBITT IAL  SQ.TJATI01TS 
9«   Seneral  Theorem.  -    The  general  homogeneous 
linear  differential  equation  is  of  the  form 


/  0 


0  .. 
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n  n-1 
d  y  d  y 

(1)  Xo   r    *  n-l      *  ^^3r  .=  0 

dx  dx 

where    Xq  ,  X^,  Xj^  are  constants  or  functions  of  x.  If 

y  =        is  a  particular  integral  of  this  equation,  then 

d  yi  d  yi 


Xq        n      +  Xj,  t       +  +    Xj^yj^  s  0 

dx  dx  - 

Multiplying  this  identity  through  hy  an  arbitrary  constant  c^, 

d  yi  cL  y^ 

(2)  XoCi  - 
dx  dx 

This  is  the  condition  that  c^yi^is  a  solution  of  ©quation  (!)• 
Similarly,  if  y  =  yg  is  another  particular  integral  of  (1), 
^-^zVz  is  also  a  solution,  and 

(3)  XoCg  +     XiCg        n-l       +  +  Xj^Cgyg  s  0 

dx  dx 

Adding  identities  (2)  and  (3) 

.  n-l 
n  n  n— 1  o.  y 

d  yi  d  <i      71   1 

XoCi   +  XoCg  5^    +  X^Ci        n-l        +  ^iCg  n-l 

dx  dx  dx  dx 


This  is  the  condition  th^^Cj_yj  +  Cgyg  is  a  solution  of 


if 

equation  (1),     Similarly,   if  n  particular  integrals  yj^ 
of  (1)  are  known. 
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(k)     y  a  c^yi  +  Cgyg  +    +  ^n^n        ^  solution. 

Since  the  general  solution  of  a  differential  equation  involves 
n  independent  arbitrary  constants,  and  since  c^,   Cg,  --- 
are  such  constants,   (k)   is  the  general  solution  of  equation  (1) 
Hence  the  theorem:-     If  y^,  yg*        -        are  n  independent 
particular  integrals  of  a  homogeneous  linear  differential 

equation  of  the  order  n,  the  function  Cj^y^  +  c^yg  +  +Cjiyn 

is  its  general  solution. 

10.   Constant.  Coefficients,     Root s  of  Auxiliary 

Equation  Heal  and  Distinct.  -     When  Xq .  X^,   X^  in 

equation  (1)  are  constants,   it  may  be  written 

d  y  d  y 

(la)  ko   ^    +  ki   ^33-    +  +  kny  »  0 

dz  dx 

Consider  the  solution  of  the  equation 

which  is  a  special  case  of  (la). 
Separating  the  variables 

+  kdx  =  0 

y 

log  y  +  kx  s  a 
log  y  -  log  c  =  -kx 


V  -kx 
JL—    =  e 

c 

-kx 

y  =  ce 


-I -J    <  f 
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This  solution  suggests  that  a  particular  integral  of  (la) 
may  he  of  the  form  y  =  e^^.     Assuming  that  this  is  the  case 
then 

,     n  mx  ^  ,      n-1  mx  _^  .  ,      mx  ^ 

k^m  e      ■♦■  ICim      e      +   +  kne      2.  0 

e°^^(kom^  +  kjm^-^  +  +  k^)  5  0 

If  m  is  a  root  of 

(6)  kom^  +  k^n^"*^  +  +  k^  *  0 

which  is  called  the  auxiliary  equation,  y  =  e'^"'^  will  he  a 

solution.     Let  m^^,  mg,  mj^  he  distinct  real  roots  of 

equation  (6)*     By  the  theorem  just  proved 

m,x  .        maX  ,  mnx 

(7)  y  =  Cj^e   ^    +  Cge  ^     +  +  c^e 

will  he  the  general  solution  of  (la) 


,2 

d  y  4y 

The  solution  of  3   ^    +  7  -  6  y  =  0    will  serve 

dx 

as  an  illustration  of  the  application  of  (7)»  Ttie  auxiliary 
equation  is 

3m*+7m-6=0 

(3m-2)(m+3)  =  0 
m  «  I      m  «  -3 
Suhstituting  in  (7) 

3*  -3x 
y  =  Cj^e      +  Cge 


1             '  ' 
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11.  Roots  of  Auxiliary  Squat 1  on  Complex  Imaginary 
all 

and  Distinct,  -     If^the  roots  of  equation  (6)  are  complex 
imaginarj'-  and  distinct,  they  will  occur  in  pairs  of  the  form 

and  u.  ~  ^i.     By  the  theorem  ^rssad  on  page  5 
solution  of  (la)  will  "be 

y  =  c^e  +      c^e  +   


+      Cn.l®     ^      A-  +       °n^  ^ 


This  can  be  written 


«^  X  /lis           ^o^x  -/lix  ,  -^.x  /fzLix  ^ 

=  c^e  ^  e  +  Cge      e  ♦  +  cn-ie  t  e  r  + 


o6„  X  -/^ix 
c  e      e  » 
n 


®  "  (^n-1®  "  °n®  ^  ) 


Jrom  Euler's  Formulas 
\  6 

e  s  Cos  <^  +  i  sin  <^ 
e        «  Cose/    -  i  sin  ^ 


V  4.*    1       /^3c  V^ix 

We  may  express  the  exponentials  e        and  e'        in  trigonometric 


form  as  follows:- 


e^  =  Cos/x  +  1  sin/x 
e  ^     =  Cos/f  X  -  i  siny^: 


t 


J  0 


■  I 


X 

1. 


Suljst ituting  these  values  in  equation  (8) 

y  ~  e^^  (cj_cos y/fiX  +  Cji  sin  +  Cgcos^^x  -  Cgi  sin^/^x) 

+  — — +  e       (c    ,  cos  /i^x  +  c     _  i  sin  /fy,x  +  c^cos 
^  n-l       /  T         n-1  /  T  ^ 

-  Cj^i  sin^i^x) 
^    l^°n-l  *  ^n^   cos^x  +  ^(^^-1  "  ^n^  sin^/^xj 


7  -  e 


<<>,x 

+  e 


If  Ci  +  Cg  =  Aj,  i(ci  -  Cg)  s  Bi,  Cg  +  C4  =  Ag,  i(c3  -  O  =  B 
and  so  on, 

(9)  y  =  ©"^^^(Ai  cos/^x  +  Bi  sin/^ix)   

+  e  ^^(A^iPos  /^x  +  B^  sin/^x) 

will  "be  the  general  solution  of  (la)  in  which,  if  Aj^,  Ag,  , 

^1*  ^2*  are  real,  c^,  C3,  must  he  assumed  to 

he  imaginary. 

Another  convenient  form  of  this  solution  may  he  obtained 
as  follows:-  Multiply  and  divide  each  term  of  equation  (9) 
hy  arbitrary  constants  kj^,  kg,  so  chosen  that 

tl  -  s±n  s.        ~  *  cos  a,,         =  sin  a^,  li.  =  cos  ag,  and  so 

—  -  Sin  a^,  X  kg  •  ^ 

on*     Such  a  selection  is  possible  for,  if  ^  =  sin  a^^,  and 

^1 

Bj  2  2  3  .Js  2 

—  «  COS  aj^,  then  A^  +  B^  =  kj^     or    kj^  -  Ik^  +  Bj^  ,  Since 
^1 

2  2 
sin  cos    a,^  *  1,  and  since 


t  ■  ■      J  r 


.0 . 


noi  J 0*3  Is 3 


either  one  of  these  numbers  may  he  taken  as  the  sin  and  the 
other  as  cos  aj^.     Similarly  for  Ag  Bg 


f     2         z~  I i 

and  so  on,     Equation  (9)   can  now  be  written 

y  =  iCje*^''^(  sin  aj^  cos/^^^x  +  cos  aj^  sinyd^x) 

k^e       (sin        cos  //gX  +  cos  ag  sin  /4x)       +  — 


2 


^/^x  +  cos  ajL  sin^/^3 


kjj^e'^^^C  sin  a^i.  cos  Z^x  +  cos  a^^  sin /^^x) 


(10)  y  =  k^e   ^     sin  {p^x  +  a^)  +  kgO        sin  r/^: 


X  +  ag)  +  


+     ky,e  ^     sin  (Z/^x  +  a^ 
/-  -J 

where  k,,  kg,  k>,  and  aj^,  ag,  .  a^,  are  arbitrary 

constants. 

To  illustrate  the  use  of  formulas  (9)  and  (10)  solve 

d  y  dy 

3   -     2         +  y  =  0 

dx 

The  auxiliary  equation  is 
3m    -  2m  +  l  =  0 


m 


2  ±'\  \  -  12  1     +  {   2  i 


Substituting  in  (9) 

y  =  e%(A-i  cos  'j  X  +  sin  j  ) 
Substituting  in  (10) 

y  s  kj_ef^sin  (t  r  +  ai) 


12,  Hoots  of  Auxiliary  Squat  ion  Multiple.  -  If 
equation  (6)  has  a  double  root  m^,  and  if  the  theorem  on  page 
holds  in  this  case,  the  general  solution  of  (la)  will  be 

nij^x            mgX  m  x 

y  =  Cj^e        +  Cge        +  c^e        +  +  c^e 

y  =  e  *  (ci  +  Cg)  +  Cge        +  +  Cne 

Ci  *  C3  is  equivalent  to  a  single  constant,  and  consequently 
only  n-1  independent  arbitrary  constants  are  involved.  This 
is  not  the  most  general  solution  for  n  such  constants  are 
necessary. 

Suppose,   that  one  of  the  roots  m^  -  m^  +  h     where  h  is 
a  small  finite  quantity.  Then 

ffiiX            (m,  +  h)x            m,x  ni_x 
^    +  c,e     ^        '     +  ce  3    +  +  cj^e  ^ 


(11)  y  s  c^e        +  Cge 


by  the  theorem  on  page   5  •     Equation  (11)  can  be  written 

(12)  y  =  e  ^  (cj  +  Cge     )  +  c^e  ^     +  +  c^^e 

hx 

Expanding  e      by  the  exponential  series  formula 


^  3    3  3  3 

e^^  -  1  +  hx  +  -f^    +  -^-^  + 


SulDstitut ing  this  value  in  equation  (12) 


3    3  3    3  -I 

y  =  e  *    Ici  +  03(1  +  hx  +  -g-y  *    ~3~y~   J 


+  Cge  +  +     c^e  ^ 


(13)     y  =  3   ^  A    +     Bx  +  Cgh  §-  X     +  cgh        X     +   ) 


mgX  m^x 
+  Cge  +  +  Cj^e 


As  the  quantities  Cj^  and  Cg  are  arbitrary  constants  we  may 
choose  *^3  ~  ^    where  B  is  any  constant.     As  h  approaches  zero 
Cgh  will  always  equal  B  and  Cg  will  "become  infinitely  large. 
Now  if  Ci  +  Cg  is  to  be  an  arbitrary  constant  A»   since  Cj^  = 
A  -  Cg,     Cj^  must  be  chosen  infinitely  large  and  negative.  Let 
h  approach  zero  in  equation  (I3)  and  substitute  A  for  Cj^  +  Cg 
and  B  for  Cgh,  Then 

(Ik)     y  =  e^i^(A  +  3x)     ^     cge^s^    +  +  c^e"^^^ 

is  the*  general  solution  of  (la). 

Suppose  next  that  three  roots  m^,  Ug,  ^3  of  equation  (5) 


are  equal.     Prom  equation  (lU)   if  = 


m 


3 


(15)     7  -  e^i^(A  f  3x)   +  cge''^^    +  +  c^e''^'' 


Let  =  ra^  +  h.     Then  substituting  in 

y  =  e^i=^(A  +  BK)  +  036^"!  *  +  +  c.e"""^^ 

jr  =  e  ^  I  (A  +  Bx)  +  CgO        +  c^e  *      +  +  Cj^e  ° 


Expanding  e      "by  the  exponential  series  formula 


^  2    2  3  3 

y  =  e^^"J7^  +  33c)  +  C3  +  cgh::  +  C3  +  C3  ^^jf-    +  __. 


+  C4'e  *       +  +  Cj^e  n 


22  33 


(16)     y  =  e'-i^f^A  +  C3)  +  (B  +  C3h)x  +  C3  +  ^3  ^ 


-J    +  c^e"'*-  +  +  c^e* 


As  in  the  preceding  case,   choose     C3  ~  where  A     is  any 

constant.     As  h  approaches  zero  Cg  will  always  equal  A3  "but 
CqIcl  and  C3  will  "become  infinitely  large.     If  B  +  Cgh  is  to  "be 
an  arbitrary  constant  Ag,    since  B  =  A3  -  C3h,  B  must  he 
chosen  infinitely  large  and  negative.     If  A  +  C3  is  to  he  an 
arbitrary  constant  A^^,   since  A  =  Aj^  -  C3,  A  must  he  chosen 
infinitely  large  a.nd  negative.     Let  h  approach  zero  in  (lo^, 
and  substitute  Aj_  for  (ii  +  C3) .  Ag  for  B  +  Cgh,  and  A3  for 

2 

C3h  ,  Then 

„    _  2  m  A  X  Tr.yi  X 

y  =  e°^i^(A,      AoX  *  A^x  )  *  c^e  *    +  >  +  CrxO 

will  be  the  general  solution  of  (la),- 


(  ,                                    "                          ■      ■  ■" 

! 

By  a  similar  line  of  reasoning  it  can  be  shown  that  if 
r  roots  of  (d)  are  equal,   the  general  solution  of   (la)  will  "be 

miX,                           2  r-1 
(17)  7  =  e  +  Agx  +  A3X    +  +  A^x       )  + 

°^r+l^  mj^x 
c^+l«   +  c^e 


Thus  the  difference  that  an  r  fold  root        produces  is 

that  the  coefficient  of  the  exponential  e        is  no  longer  an 

s  t 

arbitrary  constant,  but  a  polynomial  of  the  (r-1)  degree 
of  the  form 

2  r-l 
Jli  +  A23C  +  A3X    +  +  AyX 

involving  r  arbitrary  constants. 

As  an  illustration  of  the  application  of  (17)  consider 
the  solution  of 

3  2 

12  —5-    +  ^  —ir-   -  5  —   -2  =  0 

dx  dx  dx 


The  auxiliary  equation  is 

12  m    +Um  -5m-2=0 

1  12 
m  =  "  2   t     ni  =  —  g,     m  =  ^ 

T^om  equation  (17) 

y  =  e       (A^      Agx)  +  C3e^ 


Ik 


If  the  multiple  root  is  complex  imaginary,  its  conjugate 
will  also  "be  multiple  of  the  same  order.     Suppose  mj^  and  m^ 
of  equation  (6)  are  equal  and  complex  imaginary  of  the  form 
<^  +  /^^»  ^'^^  "^3  ^®  roots  that  are  the  conjugates 

of  m^  and  mg.     These  roots  then  will  "be  of  the  form  c?C^y/i, 
By  equation  (1?) 

y  =  e  (Ai  +  Agx)  +  e  ^  (B^  + 


Bgx) 


mgX  ™n'^ 
+  Cge        +  +  Cj,e 


This  expression  may  he  written 


y  =  e  ^  .e^    (A^  +  Agx)  +  e       .e         (Bj  +  Bgx) 


i^^s^c  m^x 
Cge        +  _  +  Cj^e  ^ 


(IS)  y  =  e   ^    [e^i    (Aj  +  A^x)  +  e  '       (B^  +  B^xW  + 


°^53f  m^x 


Cge        +  +  Cj^e  ^' 

Substituting  the  trigonometric  forms  for  &  and  e/ 

in  (18) 


o6. 
y  -  e 


"^^ l^cosy^x  +  i   sin^x)(Ai  +  Agx)  +  (cos/^x  - 

/"I              n^s^c  m  X 

ix)(3i  +  Bgx)]   +  cge        +  +  Cj^e  ^ 

(19)  y  =  ^^^f{k^  +  Bi  +  AgX  +  Bsx)   cos  /^x  +  i(Ai-Bi  +  AgX-Bgx] 


J    7           nieX  m-x 
a±n/f^xj  +  CgS^  s     +  +  c„e  ° 


Let         +  Bi  =  Ai»,         +  B3  =  Ag*,   i(Ai  -  B^)  =  A^l', 
-  Bg)   =  Ag".      Substituting  in  (I9) 
y  =  e'^^^I^Ai'  +  A2»x)cos^iX  +  (Aj"  +  Ag^x)  siny^ 


+  Cge        +  +  c  e 

n 

Trill  "be  the  general  solution  of  (la). 

Similarly  it  can  "be  shown  that  if  a  complex  imaginary 
root  is  r  fold,  the  general  solution  of  (la)  will  "be 

r-l 

Ai»  +  Ag'x  +  Ag'x"  +  +  A^»s 


(20)  y  =  e       KAi*  +  A^^x  +  Ag'x    +  +  A^'x  jcos^^x 


+  (Ai"  +  A2"x  +  Aa^x^  +  +  Aj."x^'"^)  sin/^^x^ 

"^sr+i^  ^n^ 
*  <^2r+l^   ^n« 


Another  form  of  this  result  may  "be  ohtained  as  follows:- 
writing  equation  (20)  in  its  expanded  form  and  grouping 

(21)  y  =  ©'^^''^[(Ai 

+  Ag"   sin^^ix)   +  +  x^''^(Aj.' cosy/^x  +  Aj."  sin^ 

"^sr+i^  ninX 
+  C2j.+  ie  +  +  Cj^e 

Multiply  and  divide  the  first  group  within  the  hracket  by 
kj^,   the  second  group  by  kg,  and  so  on,   the  k*s  so  chosen  that 
^i'  Ag" 

-- —    =  sin  a^,    =  sin  ag,  —    =  cos  ag,  and  similarly  for 

_=___^  ^2 


-  t 
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the  other  A* b  and  k*s.     Substituting  in  (21) 

(22)  y  *  e^^^^i  sin  {^i^      e-i)  +  ^3=^  sin  p<^x  +  ag)  +  

r-1            jl             n  nisr+i^c 
+  k^3c        sin  ^ +              "**  C2j.+  ie  +  

+  c  e 

To  illustrate  the  application  of  (20)  and  (22)  consider 
the  solution  of 

4                       3  2 

d  y           d  y           d  y  dy 

—r ^2       3+3       2+2  +1=0 

dx              dx             dx  dx 

The  auxiliary  equation  is 

4                 3  2 

m    +  2m     +  3m    +  2m  +  l  =  0 

-l+jT^i  -ll/Tl 
2                •  2 

From  equation  (20) 

y  =  e''*'^^Ai»  +  Ag'x)   cos^x  +  (A^"  +  A2"x)  sinx^cj 

From  equation  (22) 

y  «  e       Ikj^  sin  (^  x  +  a^)  +  kgX  sin  (-^  x  +  ag)  I 

HON-HOMOGENEOtrS  LINEAR  DIFFBHEFTIAL  EQUATIONS 

13«   G-eneral  Theorem.  -     The  general  form  of  the 

non-homogeneous  linear  differential  equation  is 
,n  n-1 

(1)  Xo  -4   +      -nrr   *  +  x^y  =  X 

dx  dx 

• 

1  f 

'            -  ' 

where  Xq,  X^,   -  X^  and  X  are  functions  of  x  or  constants. 

When  the  right  hand  memher  of  this  equation  is  zero,   it  has 

"been  shown  that     Y  =  y  =  Cj^y^  +  Cgyg  +  Cnyn 

general  solution,  and  for  convenience  Y  will  he  referred  to 
as  the  complementary  function  of  (1).     If  y^    is  a  particular 
integral  of  (1),  then 

n  n~l 

(2)  Xo  +  +  +  X„y^  =  X. 

dx  dx 

and  since  Y  is  the  general  solution  of  (1)  when  the  right 
hand  memher  is  zero,  then 

n  n-l 

(3)  Xo  ^  +  Xi  ^  n.i     +  +  X^Y  5  0 

dx  dx 

Adding  the  identities  (2)  and  (3) 

n                 n-1  n 
(U)  Xo  +  Xi  f^Ezi*-    *  +  X„y>   +  Xo  j4-  + 

^'~d?^  *  *  '  ^ 

This  is  the  condition  that  Y  +  y^    is  a  solution  of  (l)t  and 

-ft 

since  the  complenientary  function  contains^^ independent 
arbitrary  constants,  it  is  the  general  solution.     Hence  the 
theorem:-      The  general  integral  of  a  non-homogeneous  linear 
differential  equation  is  the  sum  of  its  complementary  function 
and  any  particular  integral, 

lU.   Symbolic  Operator  (D-m) .  -     If  the  coefficients 


Xq  ,  Xi,  in  equation  (1)  are  constants,  it  can  "be 

written 

n  n>l 

(la)  ko  f4    +        T-nrf    +  +  k^7  =  X. 

dx  dz 

,               J              s  n 
and  if  we  write     f^L  =  py;  2^        p   d_x  =  D  y, 

dz  dz^ 

equation  (la)  "becomes 

(2)   (koD""  +  kiD"""^  +  +  k^)y  =  X 

in  which  the  polynomial  ^^q^^  +  kj^D^^"^  +  .  +  kr^  is  a 

rational  integral  function  of  the  n*^  degree  "by  definition, 
and  represents  symbolically  the  differential  operator 

n  n-1 

ko  f-E-    +  ki  ^-HTT  +  +  . 

dx  dz 


m-t^^  th 


Since  any  rati onal^ function  of  the  n""  degree  with  constant 


coefficients  can  "be  written  as  the  product  of  n  linear 

factors,   some  alike,  all  alike,  or  all  distinct,   equation  (2) 

"becomes  symbolically 

ko(D-mi)(D^n2)  (D-mji)y  =  X 

By  definition  of  the  operator  D,   (D-m)y  means  -  my  and 

dz 

(D-.n)y  means     ^  -  ny.     The  product  (D-m)(D-n)y  means  (^  -  m) 

(^  -  ny)  which  equals  _  (m  +  n)  |^      mny.  Symbolically 

dz 


(3)   (D-m)(D-n)y  =  [d^  -  (m+n)  D  +  m^  y 


i 


That  is,  the  result  of  operating  on  y  with  (D-n)   first  and 
then  with  (D-m)  on  the  result  is  the  same  as  operating  on  y 
with    j^D^  -  (m+n)  D  +  Moreover,   owing  to  the  symmetry 

of  m  and  n  in  the  right  hand  memher  of  (3)»         see  that  the 
order  of  the  operators  on  the  left  is  not  essential.  Hence 
the  symbolic  representatives  of  operators  of  the  type  here 
considered  "behave  like  algebraic  quantities  for  multiplication; 
and  furthermore,  it  is  obvious  that  this  conclusion  will  hold 
regardless  of  the  number  of  operators  involved. 

15.   Operator  Method  for  Particular  Integral .  - 
Suppose  we  have  an  equation  of  the  third  order,  and  that  it 
has  been  divided  through  "by  k^.     Symbolically  it  can  be 
written 

(1)   (D-mi)(D-m3)(D-m3)y  =  X 
Let   (D-mg) (D-m3)y  =  u,  where  u  is  a  new  function. 


Then  (D-mi)u  =  X    or    |^  -  mj_u  =  X, 

This  is  a  linear  equation  of  the  first  order  and  e     ~     is  an 

integrating  factor.  Hence 

-m,x  du  V 
e     ^  -  ruiu)  =  e  X 

Integrating 

e         u  =  /  e        Xdx  +  c 
u  =  e     ye     ^  Xdx  +  ce  ^  or 

m,x 


/_       »  .         V  raiX  ^  — m^x_  m^^ 

(D-mg) (D-m3)y  =  e  ^    /  e     ^  Xdx  +  ce  ^ 


Let  (D-m3)y  =  v.  Then 

(D-m2)v  =  e  e     ^    Xdx  +  ce  ^ 

This  is  a  linear  equation  of  the  first  order  and  e"^^^  is  an 
integrating  factor.  Therefore 

'^s^  (mi-ma) 3C  /^-m,x  (t^i-t^z)^I 

ve  =  J  L®  7  ®  *  ^  dx  +  c» 

Y  ^  e    y  je  J  Q  Xdx  +  ce  Jdx  +  c»e 

/«x  m2X/'(mi-m2)x/---miX  mgX^ 

(2)v  =  e     ^e  ^e        Xdxdx    +  e     /  ce  dx  + 

c»  e 

J^^(mi-in3)x        ^         —  ^(mi-m3)x 


let  — 2   s  c".       Substituting  in  (2) 


V 


s  e  ®ye^    ^  Xdxdx    +  e      c«e     ^     ^'     +  c»e 

_  ^m2X/^(mi-ni2)x/^-miX  m  x 


nigX 

t  or 

^s^^/^  (mi-ma) x/^-miX  ^  m,x         .  ^ 

(D-m3)y  =  e    J/e'  ^  ^    Xdxdx    +  c«e   ^    +  c»e 

— mgX 

This  is  a  linear  equat  ion  of  the  first  order  and  e  is  an 

Integrating  factor,  so 

-mgx      /^f  (m2-m3)x/- (mi-m2)Xi^ -m^x  (ni-n3) 
ye  V  L®  J  ®  J/  e  Xdxdx    +  c"e 

(m,«m3)xl 
+  c'e  Idx  +  C3 

y  .  e"^3^/^^"^2-^3^V^^'"i-^2^V^-^i^Xdxdxdx+ 

^  y  m^,— mj 

+  c*  m2X  m3X 


m,x 
e  * 


21 


Let  — E   B  Ci,     —   =  Cg.  Then 

mj^-m^  1113-03 


y  -  e  ' 


Je     *        Je  Je         Xdxdxdx  +  c^e  + 


^2®        +  Cge 
is  the  general   solution  of  (1). 

th 

TTsing  the  same  argument  it  can  "be  seen  that  for  the  n 
order  equation,  the  general  solution  rill  "be 

(3)  7  •  e"'^^;^"^-!-'"''^'^/.  j;(Mx-=a)y;-»i^  X(a,)° 

+  c,e        +  Cge        +  +  c  e 

*  *  n 

in  which  the  first  line  is  the  particular  integral  and  the 
second  line  the  complementary  function. 

The  application  of  (3)   is  illustrated  "by  the  solution  of 
3  a 

3-'2*'^cLx~® 

dx  dx 

The  auxiliary  equation  is 

3  2 
m    -  3m    +  2m«0 

which  has  for  its  f irot  roots 

m^,  ~  2,  m2  ~  1|         =  0 
The  complementary  function  is 

T  =  c^e^^  +  cgs^  ^  C3 
Substituting  in  (3) 


Ox/   (l-O)xr  (2-l)xr-2x  X 
y*e/o  Je  Je      e   dxdxdx:  +  Y 


- 

•  _ 

i 

=  -yi^x  dx  +  Y 

=  _  e^(x-l)  +  Y  Pierce  -  k02 

y  =  e  (1-x)  +  cje      +  Cge    +  C3 

In  the  process  of  integration  the  constants  can  "be  ignored 
since  they  give  rise  only  to  the  terms  of  the  corapl  e^neiitary 
function. 

l6.   Variation  of  Parameters ,  -     This  method  of 
o"btaining  the  particular  integral  is  readily  applied, 
especially  if   the  order  of  the  equation  is  not  high.     It  con- 
sists in  considering  the  constants  in  the  complementary 
function  as  undetermined  functions  of  x  such  that  when  it  is 
suhstituted  in  the  given  equation  we  get  X,  and  not  zero,  as 
we  do  when  they  are  constants.     Since  we  have  n  functions  at 
our  disposal  and  only  one  condition  to  impose  upon  then,  we 
may  choose  (n-1)  other  conditions  as  we  please. 

Suppose  we  have  a  third  order  equation 

3  2 

(1)  ko  ^    +       H  +  kay  =  i 

dx  dx 

the  auxiliary  equation  of  which  has  distinct  roots.  Then 

the  complementary  function  is 

m,x  mgX  mgX 

(2)  Y  =  c^e  ^    +036        +  cge 
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t 

Differentiating  (2) 


QliZ  nisX  IQaZ  m«x  3^ 


Using  one  of  the  conditions  at  our  disposal,  let 


(3)  e»x-  dc,    ,  ac,    ,  ,-3-  003    ,  0  Then 

dz  dx  dz 


fl+)  dY  niiZ  .  mgZ  .  moX 

~  CiEij.e  ^    +  Cgmge  "    +  Cgmae 

Differentiating  (U) 

3  2  m^z  a  m  3c  2  "^3^  niiZ  a 

dx 


.  +  mge  ^  +  in3e  dc^ 

dx  dx 

Using  another  of  the  conditions  at  our  disposal,  let 

(5)  m^e"^^"^  dcj  +  mge"'^''  dcg     +  m^e'^^''  i£A    =  0.  Then 

dx  d.3c  dx 

(6)  d^Y  _  ^  „2.jniZ  2  mgx  ^  ni3Z 

Differentiating  (6) 

(7)  cl'^Y  _  ^  „  3  miZ  3  mgz  ^  3  m^x 

 Cimi  6         +  cgmg  6  ^     +  C3m3  e  ^  + 

mi  6        Zll    +  ma  e"'^^  ^^3    +  03  e  ^  ffs 


-  - 

< 

.1 . 

£10   ^  - 

1 


Substituting  (2),   (4),  (6),  and  (7)  in  (1)  we  get 

(g)  ko«ai  ®  +    ^0^2  »        -il       ^o°i3  e        *^°3     =  X 

dx  dx  4^ 

'^o'^  (3)t  ^^cL  (g)  are  three  linear  equations  in  three 

unknowns  from  which     dcj_     ,     dc^     ,  and     dc3     can  "be  found, 

dz  dx  dx 

and  c^,   C2»  and  C3  determined  "by  integration.     The  constants 

arising  from  this  integration  process  produce  the  complementary 

function  again. 

It  can  he  seen  from  this  argument  that,  if  an  n*^  order 
equation  were  givea  a  system  of  n  equations  in  n  unknowns 
would  result,  and  Cj^,  Cg,  c^^  could  he  determined. 

If  the  roots  of  the  auxiliary  equation  are  repeated  or 
complex  imaginary,  the  change  in  the  form  of  the  complementary 
function  does  not  affect  the  number  of  constants  involved, 
and  consequently  causes  no  difference  in  this  process. 

As  an  illustration  of  this  method,  let  us  solve 

(9)  H  -  3  2  ^  ' 
dx           dx  dx 

the  equation  used  to  illustrate  the  method  in  section  I5,  The 
complementary  function  is 

(10)  Y  =  c^e^^  +  Cge^  +  C3 

Differentiating 

dT      ^         2x  3C  ^     2x  X  dCa  dcg 

2JL  =  2  0,6      +  Cae     +  e  J-  +  e  — s.  +  _2. 

dx  ^  ^  dx  dx  ^ 


t 

• 

;i :;x^^vi  -  s,.  sa   

f  <?  - . 
♦ 

2.x 

Let        e      dci  +  «  '^.Cg  +  dcg     _  Then 

dx  dx 


dx 


(11)     M    =     2  cie^''  +  cg  a' 
dz 


Differentiating  this  result 


7-4    =  U  Cie23^  +  cge""  +  2  e^""  dCj^  +  e"" 

"  dx  dl" 


Let 


2x 

2  e  dc 
dx 


1    +  e 


dc 
dx 


i    a     0,  Then 


(12)  .      e^e^^  ^ 


Cge 


dx 


Differentiating  this,  we  get 


(13)  d_T 

3 

dx 


2x     ^         3C        ».     2x  ^  X  . 

g  c^e        +  Cge      +  4  e       dc^    +  e  ^<^a 

daT 


Substituting  (10),   (11),   (12),  and  (I3)   in  (9)  we  get 
2x 

(II+)     1+  e       dc^  +  e'^  = 
dx  dx 


Solving 


• 

-  ■- 

■  9  _.  ^ 
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C3  -  I  e""    +  C3 

«  Oie^''^    +     (C3  -  x)s^    +     C3  which  is  the  sane 

result  we  ©"btained  in  section  I5  "but  in  slightly  different 
form, 

17,  Undetermined  Coefficients.  -     This  nethod  for 

finding  a  particular  integral,  while  not  applicable  in  all 

cases,  is  simple  rhen  it   can  "be  used*     It  applies  to  all 

cases  in  which  the  right  hand  member  of  the  equation  contains 

only  terms  which  have  a  finite  number  of  distinct  derivatives, 

lex 

Such  terms  are     e     ,   sin  Ix,   cos  mk,  and  products  of  these, 
where  k,  1,  and  m  are  any  constants. 

The  particular  integral  is  found  by  trial.     We  are 
seeking    y  «  f(3c)  which  when  substituted  in  the  given  equation 
will  reduce  the  left  hand  member  to  the  sane  function  that 
we  have  in  the  right  hand  member.     If  y  occurs  in  the  equation, 
this  leads  us  to  choose  the  terms  of  the  right  hand  member, 
each  prefixed  by  an  undetermined  multiplier,  as  a  first  trial. 
On  substituting  this  value  in  the  given  equation  other  terms 
arise  from  the  differentiation.     Consequently,  we  select  the 
terms  of  the  right  hand  member  plus  the  terms  arising  from  its 


'   O       o    -  -J  f%  ^ 

-         ■   —        f     ■   _  .. 

-  .  -i  s  e 


Hi?  ',;o'  bo.'.- 


differentiation,  each  prefixed  "by  an  undetermined  multiplier. 
Since  this  quantity  is  to  satisfy  the  given  equation  its 
Eu"bst itut ion  in  it  results  in  an  identity.     We  may  then  equate 
the  coefficients  of  like  terms  to  determine  the  unknown 
multipliers.     This  gives  a  particular  integral  which  when 
added  to  the  complementary  function  yields  the  general 
solution  of  the  given  equation. 
Illustrative  Example, 


Solve 


«    2  +  3C 


Y  = 


c^  cos  X  +  C2  sin  X 


y>   =  ae    +  hx    -  cx  +  dx 


e 


=     ae^  +  3  hx    -  c  +  2  dx 
dx  ^ 


LSA    =  ae^  +  6  hx  +  2  d 

2 

dx 


2  ae     +  "bx    +(61)-  c)x  +  dx  + 


3  V  X  0  I 


9  - 


3.3  = 


33  :i 


L  »  P. 


2  d  -  e  =  0 
e  «  0 

SuTDstitut ing  these  values  in        and  adding  this  to  the 
complementary  function 

y  =  Cj^  cos  X  +  C2  sin  x  +  e    +  x    -  Y^: 

If  y  is  lacking,   set  X  equal  to  the  lowest  ordered 
derivative  occurring  in  the  equation  and  integrate  to  oT5tain 
y  -  f (x) ,  then  proceed  as  when  y  is  present, 

Illustrat  i  on. 

Solve  +     djr    =  X®     +     cos  X 

dX3 

-X 

Y  T  Cj  +  Cge 

=  X    +  cos  X 

dx 

3 

y  -  x_    +    sin  X    +  C3 

3 

3  2 

y^  =  ax    +  "bx    +  cx  +  d  cos  x  +  e  sin  x  +  C3 
=     3  ax    +  2'bx+c-d  sin  x+e  cos  x 


~  -^-^    =     Sax+Sh-d  cos  x  -  e  sin  x 


dx 


2 

d  yj     ^     dy^  2 

2       ♦  =     3  ax    +  (6  a  +  2  h)x  +  (2  h  +  c)  + 

dx 

(e  -  d)cos  X  -  (e  +  d)   sin  x  s  x    +  cos 


a  -  1 
3 

6  a  +  8  1)  «  0 
b  =  -1 

2  b  +  c  *  0 

c  =  2 

e  -  d  =  1 
-e  -  d  =  0 

d  =  -  1 

2 

e  =  1 
2 

Substituting  these  values  in  yy,   and  adding  to  T 

»X  3  2 

y  =  Ci  +  Cge      +1.^    -  ^    +  2x-l  cos  x  +  1_  sin  x 

3  2  2 

If  u  is  a  term  of  the  complementary  function  and  also 
of  the  right  hand  member  of 

n  ^-^ 

(1)  i^o  ^    +    1^1  — nrf    +   +  *  X. 

dx  dz 

the  method  will  fail,  for  the  substitution  of  u  or  any  of  its 

the 

derivatives  for  y  in  it  will  not  give  rise  to  u,  Prom^Theory 

t  h 

of  Squations  we  know  that  if  m  is  an  r      order  root  of  a 
given  rational  integral  equation,   then  it  is  an  ( r-1 )  order 
root  of  its  first  derivative,   (r-2)       order  root  of  its  second 

t  V) 

derivative  and  so  on  to  its  (r-l)       derivative.     If  m    is  a 
simple  root  of  the  auxiliary  equation 


1      n  ^  ,      n-1  ^  ^  ,  ^ 

kQin    +  +   ~  ^ 

it  will  not  "be  a  root  of  its  first  derivative.     This  amounts 
to   saying  that  if  u,   due  to  m,  is^parti cular  integral  of 

(2)  k,  ^    +  k,  +  +  k^y  0. 

dx  dx 

it  is  not  a  particular  integral  of  its  first  derivative,  How 
if  we  are  to  keep  u  in  the  discussion,  we  must  use  some  form 
of  it  which  on  suhstitutiom  in  (1)  will  give  rise  to  it.  This 
suggests  the  substitution  of  xu  for  y  in  (1).     This  results  in 

n                                            n-l  Ti-2 
,       d  u        ,       d  u        ,       d  u        ,    ,     ,  V   d  u 


du 

d^      ^n-1^  *  =  X 


Grouping 


4  u  ^  -._     d  u 


,n-l 

d  u  .  ,    ,     ,  V   d  u 

(^o»  T-nTT  +  ki(ji-l)       n^2     +  +  l^a  1^)  =  ^ 

dx  dx 

The  first  parenthesis  vaoiishes  "but  the   second,  which  involves 

the  derivative     of  (2) ,   does  not.     Hence  we  have  u  and  terms 

arising  from  it  "by  differentiation  and  nothing  more. 

Similarly  if  u  is  due  to  an  r  fold  root,  then  it  is  a 

particular  integral  of  (2)  and  its  first  (r-1)  derivatives, 
p 

Hence  x  u  is  the  fo  rni  we  would  select  in  place  of  u. 


To  illustrate  this  principle  for  a  simple  root  of  the 
auxiliary  equation  let  us  solve 

^    -  y  =  e 

Y  =  c^e^  +  Cge"^ 
jf,  =  Axe^  +  Be^ 

=  Axe^  +      Ae^  +  Be^ 

dx 


=  Axe^  +  2  Ae^  +  Be^ 

dx 


^    -  y,    =  2  Ae^  5  e^ 

dx 

A  =  1 
2 

X  — T         1  X 

y  =  Cj^e    +  c^e      +  2  ^ce 


For  an  illustration  of  a  multiple  root  let  ub  solve 

i!f  .  3  1^    -  2  y  =  e-^ 

dx 

-X  -X  2x 

Y  =  c,e  "    +  Caxe      +  c,e 


2    -X  -3:  _x 

s  Ax  e      +  Bxe      +  Ce 


1^  =  -Ax^e"'^  +  (2A-B)xe"^  +  (B-C)e'-'^ 


=  Ax^e-"  -  (U-3)xe-^  +  (2A-2B+C)e 

dx 


-X 


^  - 

,3  2 


-x  '   2x      1  3 

y  =  c^e  *  +  Cgxe      +  c^e      -      x  e"*^ 

This  method  will  fail  again  if  a  term  of  the  form  x^u, 
where  u  is  a  term  of  the  complementary  function,   occurs  in 
the  right  hand  memher  of  the  given  equation,     A  modification 
similar  to  the  one  used  in  the  preceding  discussion  applies 
here.     To  simplify  the  argument  suppose  that  we  have  a  third 
order  equation 

3  2 

dx  cL3c 


(1)  ko  H    +    ^1  H    *  +  kgy  =  X 


and  that  u  is  due  to  a  simple  root  of  its  auxiliary  equation, 

t +1 

The  previous  discussion  leads  us  to  substitute  x      u  for  y 
in  (1),     Substituting  and  grouping  we  get 

t+ll  ^  d  u  du  / 

(t+l)tx'*'^[^3  ^0  *  (t+l)t(t-l)x*-2/kouJ  =  X 


The  first  group  will  vanish,  "but  the  others  will  not,  and 
hence  x*u  and  terns  arising  from  it  "by  differentiation  result 
and  no  more, 

t  h 

It  can  he  readily  seen  that,   if  an  n       order  equation 

were  given,  a  similar  result  would  he  obtained. 

If  u  is  due  to  a  douhle  root  of  the  auxiliary  equation 

t +  2 

of  (1),  the  substitution  of  x      u  for  y  is  suggested  hy  the 
preceding  multiple  root  discussion.     Substituting  and 
grouping  we  get 

[3  3 
d  u  d  u  du  ~1 

d7  dx^     +    ^3  dx     +  ^3^  J  + 

V   t+lP         d%  du  "7 

(t+2)x      jj  ^0  +     2  ki  ^    +    Icguj/  + 

(t  +  2)(t+l)x*^o  1^  +  l^iuj  +  ^t+2)(t+l)tx*'^u7  =  X 

The  first  two  groups  vanish,  hut  the  others  do  not.  Hence 
t 

X  u  and  terms  arising  from  it  hy  differentiation  result,  and 
no  more, 

th 

Similarly,   if  we  have  an  n      order  equation  and  if  u  is 

1 4*7* 

due  to  an  r  fold  root,   the  substitution  of  x      u  for  y  will 
give  rise  to  x*u  and  terms  arising  from  it  by  differentiation 
and  none  other. 

As  an  illustration  let  us  solve 

dx 


f 


31  : 


no  I' 


A  -  3  X  2  . 

=  ax  e     +  Idx  e     +  cx  e 


4'  3  2 

=  (ax    +  "bx    +  cx 


dx 


^3  3  2 

=  (ax    +  bx  +  cx  +     (k  a.jc  +  3  Ijx     +2  cx)e^ 

2  -3  2  3  2 

=  (ax*  +  bx  +  cx^)e^  +  2(U  ax  +  3  bx^  +  2  cx)e^' 

dx 


2  X 
+  (12  ax    +  6  15X  +  2  c)e 


2 


-.27^  +        -  12  ax  e""  +  6  bxe""  +  2  ce     h  x  ft'' 

dx 

12  a  =  1  b  =  0 

*  '  12  ^  "  ° 

y  =  c^e     +  c^xe  e 

Summary  and  Conclusions .  -     The  problem  of  solving 
linear  differential  equations  resolves  itself  into  two  parts, 
finding  a  particular  integral  and  the  complementary  function. 

The  finding  of  the  complementary  function  when  the 
equation  has  constant  coefficients  depends  wholly  upon  the 
solution  of  an  auxiliary  equation  which  is  a  rational  integral 
function  involving  a  single  variable.     This  solution  may  beaome 
involved  if  the  order  of  the  given  equation  is  high,  but 
otherwise  no  difficulty  is  encountered.     The  form  of  the 


1 

- 

i 

r 

0. 

compleraentary  function  is  exponential  or  trigonometric 
according  as  the  roots  of  the  aii:ciliary  equation  are  real  or 
complex  imaginary. 

The  finding  of  the  particular  integral  when  the  equation 
has  constant  coefficients  may  "be  accomplished  "by  Variation  of 
Parameters,  Differential  Operator  Method,  or  Undetermined 
Coefficients,     The  first  two  methods  have  the  advantage  of 
absolute  generality,  hut,  "because  they  involve  integrations, 
may  "become  laborious  when  the  order  of  the  equation  is  high. 
•JPhe  last  method,  while  not  perfectly  general,  applies  in 
many  cases  and,   as  a  rule,   is  to  be  preferred  to  the  others 
whenever  it  can  be  used  because  it  involves  only  differenti- 
ations and  the  solution  of  simultaneous  linear  algebraic 
equations.     This  work  also  may  become  long  if  the  order  of 
the  equation  is  high,  but  usually  it  is  no  longer  than  the 
other  methods,   if  as  long. 
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